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Abstract 



In this work, we focus on the fermionic structure of the low-energy excitations of graphene (a monolayer 
of carbon atoms) to propose a new supersymmetric field-theoretic model for this physical system. In 
■ the current literature, other proposals for describing graphene physics have been contemplated at the 

level of supersymmetric quantum mechanics. Also, by observing the inhomogeneities between neighbor 
carbon atoms, Jackiw et al. have set up an interesting chiral Abelian gauge theory. We show in this 
paper that our formulation encompasses models discussed previously as sectors of an actually richer 
(supersymmetric) planar gauge model. Possible interpretations for the fields involved in the present 
graphene model are proposed and the question of supersymmetry breaking is discussed. 



1 Introduction 

Graphene is a flat, two-dimensional system, consisting of a monolayer of carbon atoms in a honeycomb array. 
Its electronic theoretical description is well known for decades [TJ [5J [3] , and its analogy to planar quantum 
electrodynamics was used in the eighties [HOIS] for the study of certain aspects such as quantum anomalies. 
Its first experimental realization, however, occurred only in 2004 [3 [8], with the follow-up experiments 
confirming the theory [5J [TO] and giving rise to a huge amount of work, a great part of which are summarized 
in various excellent introductory [TT| [TOJ [TOJ [H] or more advanced [TOJ [TOJ [T71 [TO] reviews already available. 
Thus, this (genuinely) planar carbon system appears to be a promising framework for the verification of 
ideas and methods developed in quantum (gauge-) field theories, and therefore techniques of QED 3 may 
hopefully yield a number of different and relevant results in this low-dimensional condensed-matter system 
(see reviews [TOJ [2D])- 

The fundamental electronic property of graphene [T7] is the fact that its low-energy excitations are 
described by a Dirac Hamiltonian for massless particles (called Dirac fermions). So, at this level, the only 
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field appearing in the theory is the Dirac spinor, which possesses four entries: two due to the fact that 
there are two atoms in each site of the triangular Bravais lattice (or, equivalently, there are two triangular 
sublattices performing the hexagonal one), and the remaining two due to the existence of two inequivalent 
points (the so-called 'valleys', or Dirac points, which in literature are often named K + and K_) of the first 
Brillouin zone. No considerations about the real spin of the electrons are being made here and for this reason 
such spinorial degrees of freedom in graphene carry a quantum number referred to as pseudospin. 

Another interesting issue related to graphene physics is the fact that an effective gauge field (Abelian or 
not) may also appear in the theory. This aspect is summarized in Refs. [171 119j and in the specific reviews 
[2"Tl |2"2"] . Such an effective gauge field is induced if one or more of the following possibilities are taken into 
account (besides the obvious one corresponding to an electromagnetic field) : ripples [23] or curvature effects 
of the graphene sheet; lattice distortions or strains [24] : hopping inhomogeneities between neighbour carbon 
atoms; or all of these, since one may induce the other(s). With respect to the third possibility, Jackiw 
and co-authors have EH ES [27] proposed a very interesting chiral (Abelian) gauge theory, in which also 
a complex-valued scalar field is present. The latter had been introduced previously by Hou, Chamon and 
Mudry [28] and is also generated by the hopping inhomogeneities, more precisely, by the so-called Kekule 
distortion, studied by Chamon in 2000 [29] . 

The presence of these three fields allows one to wonder if some relationship such as supersymmetry 
(SUSY) may exist among them and, eventually, other fields to be considered in addition (supersymmetric 
partners). Such a consideration is reinforced by some works [301 El 1321 1331 [Ml I3S] in which the possible 
manifestation of SUSY in a theory for graphene has already been raised, even though they remained at 
the level of supersymmetric quantum mechanics (and corresponding features of the energy spectrum), still 
lacking a proposal for a (supersymmetric) graphene field theory. 

Inspired by these results and ideas, our main contribution is to actually point out that the gauge theory 
constructed by Jackiw and collaborators EH ESI El] may be considered (at least in its original version of 
Ref. E5]) as a sector of a wider and richer planar gauge field theory: the so-called Supersymmetric r 3 -QED 
[36j . With our work, supersymmetric quantum field theory becomes one more item to be added to the list 
[19] [20] of topics originated from high-energy physics which have the potentiality to be realized in table-top 
experiments with graphene. 

The present work is organized as follows: we shall review, in the next Section, some basic facts about 
the mentioned chiral gauge theory and, in Section 3, we shall set up the main aspects of Supersymmetric 
r 3 -QED. Then, in Section 4, we shall explicitly build up a general (power-counting renormalizable) super- 
symmetric action in (2+l)-dimensional space-time to finally show that the action in Ref. |25| (and possibly 
its generalizations in Refs. [26, 27 ) could be embodied in a broader supersymmetric functional. Some con- 
clusions are depicted in the last Section. A four-subsection Appendix is added with information about the 
representation of gamma matrices. 

2 A chiral gauge theory for graphene 

In 2000, Chamon studied in detail [29] some consequences of the so-called Kekule distortion in a honeycomb 
array of carbon atoms. In 2007, Hou, Chamon and Mudry (HCM) extended this idea EE] considering a 
Kekule texture, that is, a different Kekule distortion in each point of the plane, thus introducing a scalar 
field, which will be represented here by (p. Some time later, in order to provide dynamics and finite energy 
to the vortices described by HCM, Jackiw and Pi [25] introduced a gauge field to their model, thus 
formulating the chiral gauge theory that will be summarized in the present Section. For more details about 
this theory and its extensions EHIEZ], the reader is referred to the original articles. 

HCM's original theory for the Kekule texture in graphene is described by the following Hamiltonian 
density [28] : 




where VP(r) is a four-component spinor describing the electrons in graphene (called Dirac fermions, an 
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expression meaning massless fermions) 

V i> h - ) 

where the indices a and b refer to the two triangular sublattices of the honeycomb array and the indices H 
and — refer to the two inequivalent points (Dirac points) of the first Brillouin zone. K is the 4x4 matrix 



K 



( -2id z gip{v) \ 

-2id z * g<p(r) 

g(p*{r) 2id z 

V (r) 2id z * J 



(3) 



with —2id z — j (d x — id y ) and g is a coupling constant. As mentioned above, the complex- valued scalar 
<p(v) describes the Kekule texture. Its modulus provides a single-particle mass gap for the Dirac fermions 
and its phase may describe vortices. 

Jackiw and Pi [25], concerned with the fact that HCM's model leaves unspecified the dynamics that gives 
rise to the complex vortex profile, proposed an extension of HCM's model introducing a gauge potential and 
coupling it to the Dirac fermions in a chiral manner. 

Using the Dirac matrices forms 



a = (a 1 ,a 2 ,a 3 ) = (° _°J 
and the chiral gamma matrix, as the Hermitian one 
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the gamma matrices are 



1j = Pol 



1°j 
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With these matrices, K in Eq. ([3]) may be represented as 

$ t M = $ t («.p + g/3 [ip r - iip* 7 5 J ]) 



(4) 



(5) 



(G) 



(7) 



where p is the operator — iV and tp = ip r + . Thus, the final Dirac Hamiltonian density, with the additional 
gauge potential, can be written as 



* f K A * = *t a ■ [p - q~f( A] * + g^/3 [ip r - i 7/ ip*]^ 

= 7 j • (p - qA) + 7 j ■ (p + qA)f_ + g ip *+ + gip* * H 



(8) 



where q is a coupling constant describing the chiral charge (in Ref. [25 , q is set to unity), the Dirac adjoint 
is ^ = ^ 7° and the chiral components are ^± = | (1 ± 7/) The Hamiltonian density which is present 
in Eq. J51) is invariant under a local chiral transformation 



e 2lqu cp 



A-> A + Vw 



(9) 
(10) 



Notice that this system possesses a global fermion number symmetry, with just the Fermi fields trans- 



forming with a constant phase: ^ 



\Jj. Consequently, the theory possesses a local chiral U(l) symmetry 
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and a global U(l) fermion number symmetry 0. Because the theory resides in (2 + 1) dimensions, no chiral 
anomalies interfere with the chiral gauge symmetry. 

With the additional gauge potential A, the Dirac eigenvalue problem in this model differs from HCM's. 
According to Eq. ® , 

[a ■ (p - qi ( A) + g(3(<p r - 1 7/ p*)]* = E^l . (11) 

Notice that a 3 , which will now be renamed as R, anti-commutes with the matrix structure on the left side of 
Eq. (fTTj) . Therefore if ^e is an eigenfunction with eigenvalue E, R^e belongs to eigenvalue — E, and zero 
modes can be chosen as eigenstates of R. This is a consequence of the "sublattice symmetry" identified by 
HCM and renamed as energy-reflection symmetry in Ref. [27] . 
From Eq. (|8[), we can write down the Lagrangian density 

£hcm-jp = *+7$(vD+)*+ + *-7j(^)*- - .9^*+*- - .9^**-*+ , (12) 

where iD^ — id^ — qA^ and iD~ — id^ + qA^. 

To set up this Jackiw-Pi's graphene theory in the context of supersymmetric T3-QED, which will be 
presented in the next Section, it is convenient to re- write the Lagrangian density, given in Eq. (|12p , in terms 
of two-component spinor notation. The four-component spinors ^, ^+ and can be written in terms of 

f tp b +\ ( V-\ 

two-component spinors ip+ — and -0- = 1 80 that 



*= : . *+= : ), *-= „ ■ as) 





We wish to keep the symbol ip- (without "prime" ) for future spinor redefinition. The 4x4 Dirac matrices 
can be expressed in terms of the minimal 2x2 Dirac matrices as 

u _ ( 7^7° 



^ = { W ) (14) 

with n spanning the values (0,1,2) and 7 M = (a z ,ia y , —i<J x ). In the sequel, using Eq. (jlU)) and Eq. (IT4")) . we 
can write down £hcm-jp in terms of two-component spinors and minimal 2x2 Dirac matrices as 

£hcm-jp = tp + ^ IJ, {iD^)ip + + il>_^{iD^)ip- - gpip + ip- - gtp*ip_ip+ , (15) 

where 

V>-=7V1 = I I • (Ki) 




3 The Supersymmetric T3-QED Action 



With the purpose of formulating the supersymmetric T3-QED action, we refer to the work by Salam and 
Strathdee [37], where the superspace and the superfields in (3+1) space-time dimensions were originally 
introduced. Extending these ideas to the present case, the elements of the superfield are parameterized 
as (x^,0), where are the coordinates of the space-time and the fermionic coordinates, 9, are Majorana 
spinors, 9 C ~9. Now, we are ready to introduce the (simple) supersymmetric t 3 -QED formulation by means 
of the superfield formalism. As a first step, we define the complex scalar N=l, D=(2+l)-supcrfields with 
opposed [/(l)-charges, <!>+ and as 

$± = A ± + 9^ ± - -66F± and & ± = A* ± + ~$ ± 6 - ^99F* ± , (17) 

1 The gauged Abelian symmetry does not represent the usual electric charge — gauge boson interaction. As a matter of fact, 
it is the additional global Abelian phase symmetry that happens to be related to the electric charge. 
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where A± are complex scalars, ip± are Dirac spinors and F± are auxiliary scalar fields that obey the following 
supersymmetric transformations: 

5A ± = eV± (18) 
5ip± = sF± + ie^d^A± (19) 
5F± = iTfd^ . (20) 

Notice that, by introducing A± and ip± together in the same superfield, the scalars cannot be neutral 
under the global fermion number U(l) symmetry [25 | I27 j mentioned in the previous Section. 

In (2+l)-D, N=1-SUSY accommodates neutral matter in real scalar supcrficlds. To describe charged 
fields, we have to combine two real supermultiplets. 

In the Wess-Zumino gauge, the gauge superconnection, r a , is written as 

Ta = i (7^) a A v + SOK and r a = -i (9^) a A^ + 96\ a , (21) 

where A^ is the gauge boson and A a is its partner, the gaugino (Majorana spinor). Defining the "field- 
strength" superfield W a as 

W a = -~D b D a T b , (22) 

with covariant derivatives given by 

D a = d a - i (r0) a and D a = -d a -i (9^) a 0„ , (23) 

we obtain 

W a = A Q + TfthFp, - l -e9^ ab (d^X b ) (24) 

and 

W a = A a - AKbF.u - l -ee (a„X 6 ) 7 £ 6 , (25) 

where = \ [7^,7"] are the generators of the Lorentz group in (2+1) space-time dimensions. 

The gauge covariant derivatives that act on the matter fields with opposed t/(l)-charges, $+ and 
are respectively given by 

V$± = (D a T iqT a ) $± and = (A, ± iqT a ) $L • ( 2 6) 

Using the definitions previously given for the superfields in Eqs. (fTT)) . (I2"TT) . ([M]), (|2"5j) and the gauge 
covariant derivatives given in Eq. (|26l) . we can build the supersymmetric T3-QED action [36]: 

5 r3 _ Q ED = J d 3 xd 2 6 l.~WW + (V$^)(V$+) + (V«l) (V*_) + 2m - | . (27) 

In terms of component fields, the action given by Eq. (|27l) may be written as 

SV 3 -QED = J cPx^iJ-fd M X-^F MV F' w + 

- A* + UA+ - A* QA_ + iij + ^d^ + + Up^d^- + F^F+ + F1F- + 

- qAf, (i> + ^iP + - ^p_Y l ^P- + iA* + d^A+ - iA*_d^A_ - iA+d^Al + iA-d»A*_) + 



- iql A + tp + X - A^tp_X - A* + Xip+ + + q A^M 1 \A\A+ + A*_A_ j + 

-m(i} + tp + -i>_ip- + A* + F+ -A*_F_ + A+FX -A_F* ] 1 . (28) 



As one can easily read from Eq. (|28p , the two-flavoured fermionic sector presents diagonal minimal 
coupling with the gauge boson, in correspondence with the first two terms displayed in Eq. (|15[) . Nevertheless, 
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the Yukawa-like interaction sector of Eq. (|15|) does not show up in the above T3-supersymmetric Lagrangian, 
even if we take the auxiliary fields F + and F_ on-shell. As a matter of fact, a 3> 4 -sector has to be supplemented 
in order to provide the theory with the proper vertices, and such an extension must adopt the symmetry 
content of Lagrangian (f28f as a guideline for its construction. Thus, we re-state the invariance of the 
supersymmetric T3-Lagrangian with respect to an Abelian gauge transformation associated to the potential 
and to the field-strength F^ u , and we render evident a symmetry with respect to the following discrete 
parity variation: 

$± -> -$ T , (29) 

as space-time transforms according to x 1 , x 2 ) — > (x°, —x 1 ,! 2 ). 

The superfield variation expressed in Eq. (|29|) has its component-field counterpart spanned as 



i>' ± = tViM (30) 
A± = -A T ; 

F'± = Ft, 



where the 4-component spinor representation of parity reads 



with r Par defined in the context of a 4-spinor Pisarski representation, related to the Jackiw representation 
as drawn in the Appendix. 

To complete the variation rule for the supersymmetric action, we shall comment on the fact that the 
Grassmannian (2+l)D-integration measure has an odd character with respect to parity, 

d 2 e = dSde ^> d 2 e' = - d 2 e. 

This stems from imposing homogeneity on the effect of parity action onto superfields $±. The scalar 
superfields <$>± host a contracted two-component fermionic term in their expansions, namely, 6ip±. As ift± 
transforms according to Eq. (f5T7|) . 9 is required to change tc^| 7°7 2 9, up to a phase, as the term (6ip±)' must 
still be a scalar. So, regardless of the phase, d 2 9 is necessarily odd under parity. As a consequence, the 
superpotential to be added to the theory has to be odd as well (note that the bare mass term proposed in 
Eq. (|27|) already shares such a property). 

4 The $ 4 -sector 

The most general local U{\) and parity-invariant supersymmetric action with a $ 4 -sector reads: 

= jjd 3 xd 2 9 ($^$+) 2 - (V_$_y 



+g J d 3 xd 2 9 

Sf + Sg, 



(31) 



where / and g are (real) coupling constants (to be associated later with Jackiw-Pi's g present in Eq. (fTS]) ). 
The quartic term with coupling constant g explicitly breaks the global U(l) fermion number symmetry. 
Nevertheless, we decide to keep this term in order to describe possible regimes where the chiral interaction 
dominates over the forces dictated by the global symmetry. The chiral symmetry has a dynamical character 



2 Such a variation also matches consistency as one considers both super-translations, — > + ie-y^O; 9 — > 8 + e, and parity, 
x' 1 — > —x 1 ; x°,x 2 left unchanged. Also, the one-flavor content of the Grassmannian sector of superspace, as one deals with an 
N=l theory, requires the variation of 9 to be associated to a 2x2 block of r Par . As a room for further extensions, we comment 
on the fact that the model here proposed is obtained after a truncation of an N=2,D=(2+1) descent of a N=l,D=(2+2) 
supersymmetric theory 36 . 
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in that it dictates a gauge interaction; on the other hand, the global fermion number appears to have only a 
kinematic character. The latter classifies the states and held configurations (e.g. vortices in Refs. [25j [27]) 
without, however, introducing gauge-type interactions. Nothing prevents us from setting g = whenever we 
wish to recover results for which fermion number conservation is mandatory. 
The component-wise $ 4 -action expression follows: 



s t 



- 1 d 3 xC f , 



with 



{F + A* + +f;A + ) \A + \ 2 - (F_A*_+F1A_) L4_| 2 + 2(L4 + |V + ^ + - |A_| 2 ^_) + 



+ - (Ali> + i> c + + Afr+i>+ ~ A 2 _^^ c _ - A* 2 ^-) 



(32) 



Also, S g = J d 3 xCg hosts the Lagrangian density 



A+\ z [ A_F+ + A*_Fl + ^( A + F - + A *+ F -) + V4^- + ^-V>+ ) + 



- \A-\ A yA + F- + A* + F1 + ^{A-F+ + A*_F+) + + ^+V>- ) ■+ 

+ \ {J$$ + il>- + A* 2 ^_^+ - A 2 Jp_^+ - A*_ 2 V^-) + 
+ (A- A+ + A*_A* + ) (^+ - + 
+ i [A* + A- - «-) + ALA+ ~ + 

+ i (A 2 + A-F% + AfA*_F+ - A 2 _A+F*_ - A* 2 A* + F-) , 



(33) 



where the two-component charge- conjugated spinors are defined by ip^_ = ia x ip+ and = ia x ip*_. 

As we aim to end up with an action on-shell for the auxiliary sector, the field equations for F± are 
obtained: 



SS 
SF+ 


= F+ 


— mA*j_ 


-fA* + \A + \ 2 


~9 


\A + 


| 2 A_- 


\- \ (A* 2 A*_ 


- |>4- 


\ 2 AJ) 


= 


SS 
SFZ 


= F_ 


+ mA_ 


+ fA_\A_\ 2 


+ 9 


\A- 


?A\- 


t- \ (A 2 -A + 


-|A+ 


\ 2 A* + ) 


= 


SS 


- F+ 


— mA+ 


-fA+\A + \ 2 


- 9 


\A + 


\ 2 A*_ - 




- 1-4- 


\ 2 A*_) 


= 


SS 

SF- 


= Fl 


+ mA*_ 


+ fA*_\A-\ 2 


+ 9 


\A- 


\ 2 A + - 




-1^ 


-\ 2 A + ) 


= 



(34) 

The outcome renders the physical scalar sector of the complete action (that results from summing up 
Eq. ([IT]) and Eq. $T$) as follows: 



V = F^F + + F1F_ = 
= m 2 (\A+\ 2 + \A\ 



2m/(L4 + | 4 + |A_| 4 



rag 

1 2 1 . ,2 



(\A+\ 2 + |A_| 2 ) (A+A- + A* + A*_ 



\A-\*j {A+A- + A* + A*_) - fg\A+\ 2 \A-\ 2 (A+A_ + A* + A* 



AJ 



-^L4 + | 2 L4_| 2 ( I 



4 



\A+\ Z + |A_H \A\A 2 _ + A\ 2 A*_ 2 



I 2 + \A-\ 2 



»-[\A+f + \A " 



(35) 
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where a minus sign has been omitted to allow for the direct identification of the potential V. 

At this point, we would like to reassess the parity-symmetric guidance we followed to propose the $ 4 - 
sector. As a matter of fact, the mass parameters present in extended models that follow original Jackiw-Pi's 
chiral gauge theory [26] [27] stand for v.e.v.'s of scalar fields, and a particular field is taken to be odd with 
respect to parity. Coherently, the corresponding mass parameter flips its sign as parity acts upon fields. 
If we had assumed the parameter m displayed in Eq. (|2T| to be odd as well, the parity transformation 
here proposed would not leave the action invariant, as it would also change by a sign. Instead, we take 
the viewpoint of considering the mass parameter as a fixed "bare" mass, here included for the sake of 
generality, left unchanged by parity as the physical mass would get contributions from non-trivial minima of 
the power-counting renormalisable A 6 -potential (mind Eq. ((35])). The corresponding critical values of proper 
combinations of A-fields should inherit the odd-parity property proposed by Jackiw et al. Also, as far as the 
parameters / and g are concerned, those are taken invariant under parity as well, for the premise of their 
identification with v.e.v's of scalar fields would lead to a primitive non-renormalisable ^-potential, which 
we decided to avoid. 

For the simple case g = 0, the associated potential bears four minima: 



(0, 0), trivial case (36) 




where the entries refer to values of \A + \ and |j4_|, respectively. 

If we consider the mass matrix for the whole set of original fermionic fields (A, ip+, tjj—, ipt.) and 
evaluate its eigenvalues on the asymmetric (|38[) minimum, we are left with three degenerate null outcomes 
and two degenerate — m results. So, the simplest self-interaction scenario already offers the possibility of 
describing massive Dirac fermions in graphene, a fact associated to textures and other condensed matter 
systems. Massive charged fermions in graphene, after supersymmetry breaking takes place through some 
mechanism (we do not address to this issue here), might combine into bound states of a scalar character and 
their attractive interaction may reveal non-trivial and interesting effects of the doubly-charged scalars. 



5 Conclusion and Perspectives 



Back to the component-wise expression for the spinor-scalar 3- and 4-vertex sectors of the complete super- 
symmetric action, we write 



afx { -iql A + ip,X - A^ip_X - Al \ip+ + A*_Xtjj. 



2/IA, 



g {A + A_ + A* + A*_) 



2/|A_| 2 + g (A+A- + A* + A* 



1p_7p— — 



! -A\ 2 + 9 -A\A 



f_ A 2 , 9 



—A* 2 + -A* A 
2 ~ 2 ~ 



f A2 , 9 



—A 

2 



2 + 



2 + 2 



A+A* 



+ g\A.\ 2 (^++^ + r-) -§(A^-Ai 2 ) </>+V- -f {a% 2 -A 2 _)^_^ + ) 



(40) 



The content of the above Lagrangian density together with the fermionic minimal couplings with the 
gauge boson displayed in Eq. ([28]) demonstrate that the supersymmetric T3-QED theory given in Eq. (f27]) 
supplemented by a <I> 4 -term given in Eq. (|31j) provides a theoretical framework that extends Jackiw-Pi's 
original chiral gauge theory [3S] (represented here by Eq. (TT5]) ). Furthermore, room is also available to 
allow for accommodating extended models 26, 27] in a supersymmetric scenario, with the expectation that 
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it may describe some of the physical features of graphene. Of course, some identifications are necessary: 
for example, in a crude comparison, one would identify Jackiw-Pi's scalar field ip with squares of the type 
^4+, A 2 _ (and complex conjugates) which arise from our quartic superfield action, Eq. (pTTj) . However, to 
exactly identify terms and degrees of freedom one should establish the proper combinations of fermions 
that diagonalize the mass matrix (eigenvectors) upon a particular choice of scalar fields configuration that 
minimizes the potential. The study of critical points for the complete potential (Eq. ([35]) ) as well as the 
analysis of mass matrix eigenvectors will soon be reported elsewhere. 

We would like to stress — and this is a non-trivial remark — that the supersymmetric T3-QED related to 
the chiral gauge theory is nothing but a planar descent of an (N=l)-supersymmetric Abelian gauge theory in 
(2+2) space-time dimensions [36]. It remains to be understood — and, hopefully, we shall also be clarifying 
this issue soon — how to relate graphene physics to the dimensional reduction of a gauge theory set up in 
Atiyah-Ward space. This may be a rich problem to be understood. 

A comment on the gaugino is in order. It is an additional field with respect to the work of Ref. [25) . 
To properly interpret the physical role of this (neutral) fermionic field, one has to consider that it mixes 
with the charged fermions whenever the complex scalars acquire non-trivial vacuum expectation values, as 
shown at the end of the previous Section. In this case, the A-field itself does not describe a mass eigenstate 
and, therefore, it cannot be interpreted as a physical excitation. If we choose a particular vacuum, one 
has to diagonalize the mass matrix in order to identify the correct combinations of the fermionic fields that 
correspond to the mass eigenvalues. Actually, one has to check if SUSY is also broken and identify which 
combination of fermions describe the Goldstino, the massless neutral excitation associated to the spontaneous 
breaking of SUSY. To the best of our knowledge, these properties cannot yet be accommodated into the 
phenomenology of graphene. The physics of graphene has a very rich phenomenology. Supersymmetry may 
indeed open up a number of interesting issues in connection with the couplings between the fermionic and 
bosonic excitations that may arise in connection with this monolayer system. The spontaneous breaking 
of N=l-supcrsymmetry in (2+l)-D is less trivial than the same question in (3-(-l)-D, since N=l does not 
exhibit a complex structure underneath, as planar N=2-supersymmetry does. So, we conclude this work with 
the perspective to reassess the problem of spontaneous supersymmetry breaking in (2+l)-D by adopting 
graphene as a physical system to apply our results. In particular, supersymmetry breaking should help us 
to understand and classify the rich structure of couplings in graphene. 
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A Appendix 

A.l General Considerations 

All the gamma matrices respect the anticommutation and Hcrmitization relations respectively given by 



where rf v = diag(l, — 1, —1). The other gamma matrices will always have the following properties: 



(A.l) 



75 2 = 1, 75 f = 75 ; 



7' 



1, 7 3T = "7 3 



(A.2) 



and the additional anticommutation relations given by 



{l* 1 , 75} = , 



ir , 7 3 } = , 



{7 3 , 75} = . 



(A.3) 
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A. 2 Jackiw-Pi representation of gamma matrices 

The 4x4 gamma matrices in Jackiw-Pi rep. [25] . which will be called 7j, can be expressed in terms of 
gamma matrices in minimal 2x2 representation, which we will represent by 7 M , as 

u_f 7^7° 



* = { W ) (A ' 4) 

with /x spanning only the values (0,1,2) and 

7 M = {<J z ,ia y ,-ia x ) . (A. 5) 

The additional matrices 7/ and 7}, are given by 

I 1 °\ .,3 _„-„,0„l„2„,J_ / 



7 5 J =^o 75 = = ^ > 'oj- (A.6) 

Since 7j is not associated to the space-time, it will behave as a second 7 5 -type matrix, that is, it anticommutes 
with all the gamma matrices associated to the space-time. 

A. 3 Pisarski representation of gamma matrices 

The 4x4 gamma matrices in Pisarski rep. [38], which will be called jp, can be expressed in terms of gamma 
matrices in minimal 2x2 representation, which we will represent by 7^, as 



-7" 

with \i spanning only the values (0,1,2) and 

7 M = (a z ,ia y ,-ia x ) . (A 
The additional matrices 7^ and jp, are given by 



>:*=( °_j 0), a. = Hhhf^ = ( ? )■ (A ' M 

Since jp is not associated to the space-time, it will behave as a second 7 5 -type matrix, that is, it anticommutes 
with all the gamma matrices associated to the space-time. 

A. 4 The transformations that relate Jackiw-Pi. rep. to Pisarski rep. 

With the aid of the matrix S = ( _ n ], we can write 

7 / 

(A.10) 
(ATI) 
(A.12) 

Eq. (|A. 12|) follows from jp and jj definitions given, respectively, in Eqs. (|A.9[) and (IA.6I) and, also, 
from the transformations given in Eqs. (jA.101) and (| A. 1 1|) . In the sequel we will write down the inverse 
transformations, 

1'j = S lp ^S (A. 13) 

7/ = Sjp^S (A. 14) 

7,5 = S 1P S . (A.15) 



7^ 




7 5 P 




7^ 


= S 7 3 jS. 
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In this work, we focus on the fermionic structure of the low-energy excitations of graphene (a monolayer 



Abstract 



of carbon atoms) to propose a new supersymmetric field-theoretic model for this physical system. In 
the current literature, other proposals for describing graphene physics have been contemplated at the 
level of supersymmetric quantum mechanics. Also, by observing the inhomogeneities between neighbor 
carbon atoms, Jackiw et al. have set up an interesting chiral Abelian gauge theory. We show in this 
paper that our formulation encompasses models discussed previously as sectors of an actually richer 
(supersymmetric) planar gauge model. Possible interpretations for the fields involved in the present 
graphene model are proposed and the question of supersymmetry breaking is discussed. 

• rH 

X 

1 Introduction 

Graphene is a flat, two-dimensional system, consisting of a monolayer of carbon atoms in a honeycomb 
array. Its electronic theoretical description is well known for decades [H [21 12] , and its analogy to planar 
quantum electrodynamics was used in the eighties [H [5J [5] for the study of certain aspects such as quantum 
anomalies. Its first experimental realization, however, occurred only in 2004 [7] [5], with the follow-up 
experiments confirming the theory [UJ [TU] and giving rise to a huge amount of work, a great part of which 
are summarized in various excellent introductory [TTJ [12] |T3j [M] [15] or more advanced [TU [TTJ [TH] [19] [20] 
|2"T1 l2"2l |2"51 |2"T1 125] reviews already available. Thus, this (genuinely) planar carbon system appears to be a 
promising framework for the verification of ideas and methods developed in quantum (gauge-) field theories, 
and therefore techniques of QED3 may hopefully yield a number of different and relevant results in this 
low-dimensional condensed-matter system (see reviews [2TH 127] ). 



* E-mail: evertonabreu@ufrrj.br 
'E-mail: marco@cbpf.br 
-^E-mail: lpgassis@ufrrj.br 
§E-mail: helayel@cbpf.br 
^E-mail: nogue@cbpf.br 
II E-mail: paschoal@cbpf.br 
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The fundamental electronic property of graphene [20] is the fact that its low-energy excitations are 
described by a Dirac Hamiltonian for massless particles (called Dirac fermions). So, at this level, the only 
field appearing in the theory is the Dirac spinor, which possesses four entries: two due to the fact that 
there are two atoms in each site of the triangular Bravais lattice (or, equivalent ly, there are two triangular 
sublattices performing the hexagonal one), and the remaining two due to the existence of two inequivalent 
points (the so-called 'valleys', or Dirac points, which in literature are often named K + and K_) of the first 
Brillouin zone. No considerations about the real spin of the electrons are being made here and for this reason 
such spinorial degrees of freedom in graphene carry a quantum number referred to as pseudospin 0. 

Another interesting issue related to graphene physics is the fact that an effective gauge field (Abelian 
or not) may also appear in the theory. This aspect is summarized in Refs. [20] [26] and in the specific 
reviews [301 EI] ■ Such an effective gauge field is induced if one or more of the following possibilities are taken 
into account (besides the obvious one corresponding to an electromagnetic field): ripples |32] or curvature 
effects of the graphene sheet; lattice distortions or strains — see, e.g., Ref. [33]; hopping inhomogeneities 
between neighbour carbon atoms; or all of these, since one may induce the other(s). With respect to the 
third possibility, Jackiw and co-authors have [34J [35] [36] proposed a very interesting chiral (Abelian) gauge 
theory, in which also a complex- valued scalar field is present. The latter had been introduced previously by 
Hou, Chamon and Mudry [37j and is also generated by the hopping inhomogeneities, more precisely, by the 
so-called Kekule distortion, studied by Chamon in 2000 [38] . 

The presence of these three fields allows one to wonder if some relationship such as supersymmetry 
(SUSY) may exist among them and, eventually, other fields to be considered in addition (supersymmetric 
partners). Such a possibility is reinforced by some works [39] [40] [41] [42] [43] [44] in which the possible 
manifestation of SUSY in a theory for graphene has already been raised, even though they remained at 
the level of supersymmetric quantum mechanics (and corresponding features of the energy spectrum), still 
lacking a proposal for a (supersymmetric) graphene field theory. 

Inspired by these results and ideas, our main contribution is to actually point out that the gauge theory 
constructed by Jackiw and collaborators (34] [35] [36] may be considered (at least in its original version of 
Ref. [34]) as a sector of a wider and richer planar gauge field theory: the so-called Supersymmetric t 3 -QED 
[45] . With our work, supersymmetric quantum field theory becomes one more item to be added to the list 
[26] [27] of topics originated from high-energy physics which have the potentiality to be realized in table-top 
experiments with graphene. 

The present work is organized as follows: we shall review, in the next Section, some basic facts about 
the mentioned chiral gauge theory and, in Section 3, we shall set up the main aspects of Supersymmetric 
T3-QED. Then, in Section 4, we shall explicitly build up a general (power-counting renormalizable) super- 
symmetric action in (2+l)-dimensional space-time to finally show that the action in Ref. [34] (and possibly 
its generalizations in Refs. [35, 36 ) could be embodied in a broader supersymmetric functional. Some con- 
clusions are depicted in the last Section. A four-subsection Appendix is added with information about the 
representation of gamma matrices. 

2 A chiral gauge theory for graphene 

In 2000, Chamon studied in detail [35] some consequences of the so-called Kekule distortion in a honeycomb 
array of carbon atoms. In 2007, Hou, Chamon and Mudry (HCM) extended this idea [37] considering a 
Kekule texture, that is, a different Kekule distortion in each point of the plane, thus introducing a scalar 
field, which will be represented here by (p. Some time later, in order to provide dynamics and finite energy 
to the vortices described by HCM, Jackiw and Pi [34] introduced a gauge field to their model, thus 
formulating the chiral gauge theory that will be summarized in the present Section. For more details about 
this theory and its extensions [36] , the reader is referred to the original articles. 

HCM's original theory for the Kekule texture in graphene is described by the following Hamiltonian 
density [37] : 




(1) 



1 Nevertheless, very recent works 28 29 present strong evidence (both experimental and theoretical) that pseudospin must 
be interpreted as a real angular momentum. 
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where ^(r) is a four-component spinor describing the electrons in graphene (called Dirac fermions, an 
expression meaning massless fermions) 

v ^ / 

where the indices a and b refer to the two triangular sublattices of the honeycomb array and the indices H 
and — refer to the two inequivalent points (Dirac points) of the first Brillouin zone. K is the 4x4 matrix 



(2) 



K = 



( -2id z gip( r ) \ 
-2id z * g<p(r) 

gip*(r) 2id z 



V g<p*(r) 2id z 



(3) 



/ 



with —2id z = i (d x — id y ) and g is a coupling constant. As mentioned above, the complex- valued scalar 
f(r) describes the Kekule texture. Its modulus provides a single- particle mass gap for the Dirac fermions 
and its phase may describe vortices. 

Jackiw and Pi [33], concerned with the fact that HCM's model leaves unspecified the dynamics that gives 
rise to the complex vortex profile, proposed an extension of HCM's model introducing a gauge potential and 
coupling it to the Dirac fermions in a chiral manner. 

Using the Dirac matrices forms 



a = (a\a 2 ,a 3 ) = Q 
and the chiral gamma matrix as the Hcrmitian one 



P = 



/ 
/ 



it 



12 3 

-ia a a 



(7 5 J ) 2 = I, 



the gamma matrices are 



1j = Pol 



7°=/3, 7 5 J = ^°jlhhj- 



With these matrices, K in Eq. ((3j may be represented as 



* f K * = * f (a ■ p + gf3 [tp r - iip 1 7/]) 



(4) 



(5) 



(6) 



(7) 



where p is the operator — iV and ip 
gauge potential, can be written as 



ip r + i(p l . Thus, the final Dirac Hamiltonian density, with the additional 



* f K A * = * f a ■ [p - q~f( A] * + g^/3 [ip r - i 7/ tp^ty 

= 7j • (p - ffA) *+ + 7.7 ' (P + gA)tf _ + g ip 



(8) 



where q is a coupling constant describing the chiral charge (in Ref. [34], q is set to unity), the Dirac adjoint 
is ^ = W 7° and the chiral components are ^± = | (1 ± 7/) \&. The Hamiltonian density which is present 
in Eq. (|5J) is invariant under a local chiral transformation 



e 2 ^ 95 



$ _> e *9"7 5 $ , A A 



Vu ; 



(9) 
(10) 



Notice that this system possesses a global fermion number symmetry, with just the Fermi fields trans- 



forming with a constant phase: '3/ 



VP. Consequently, the theory possesses a local chiral U(l) symmetry 



A Super symmetry model for graphene 



4 



and a global U(l) fermion number symmetry 0. Because the theory resides in (2 + 1) dimensions, no chiral 
anomalies interfere with the chiral gauge symmetry 

With the additional gauge potential A, the Dirac eigenvalue problem in this model differs from HCM's. 
According to Eq. ® , 

[a ■ (p - qi ( A) + g(3(<p r - i 7/ ^)]* = EV . (11) 

Notice that a 3 , which will now be renamed as R, anti-commutes with the matrix structure on the left side 
of Eq. (fTT|) . Therefore, if is an eigenfunction with eigenvalue E, R^>e belongs to eigenvalue —E, and 
zero modes can be chosen as eigenstates of R. This is a consequence of the "sublattice symmetry" identified 
by HCM and renamed as energy-reflection symmetry in Ref. [3 6) . 
From Eq. (|8|), we can write down the Lagrangian density 

£hcm-jp = + *-7j(^)*- - .9^*+*- - <7^**-*+ , (12) 

where iD^ — id^ — qA^ and iD~ — id^ + qA^. 

To set up this Jackiw-Pi's graphene theory in the context of supersymmetric T3-QED, which will be 
presented in the next Section, it is convenient to re- write the Lagrangian density, given in Eq. (|12p , in terms 
of two-component spinor notation. The four-component spinors ^, ^+ and can be written in terms of 

f tp b +\ ( V-\ 

two-component spinors ip+ — and = , so that 



*= : . *+= : ), *-= „ ■ as) 





We wish to keep the symbol ip- (without "prime" ) for future spinor redefinition. The 4x4 Dirac matrices 
can be expressed in terms of the minimal 2x2 Dirac matrices as 

u _ ( 7^7° 



^ = { W ) (14) 

with n spanning the values (0,1,2) and 7 M = (a z ,ia y , —i<J x ). In the sequel, using Eq. (jTB")) and Eq. (IT4")) . we 
can write down £hcm-jp in terms of two-component spinors and minimal 2x2 Dirac matrices as 

£hcm-jp = tp + ^ IJ, {iD^)ip + + il>_^{iD^)ip- - gpip + ip- - gtp*ip_ip+ , (15) 

where 

V>-=7V1 = I I • (Ki) 




3 The Supersymmetric T3-QED Action 



With the purpose of formulating the supersymmetric T3-QED action, we refer to the work by Salam and 
Strathdee [46] . where the superspace and the superfields in (3+1) space-time dimensions were originally 
introduced. Extending these ideas to the present case, the elements of the superfield are parameterized 
as (x^,0), where are the coordinates of the space-time and the fermionic coordinates, 9, are Majorana 
spinors, 9 C ~9. Now, we are ready to introduce the (simple) supersymmetric t 3 -QED formulation by means 
of the superfield formalism. As a first step, we define the complex scalar N=l, D=(2+l)-supcrfields with 
opposed [/(l)-charges, <!>+ and as 

$± = A ± + 9^ ± - -00 F± and & ± = A* ± + ~$ ± 9 - \^99F* ± , (17) 



2 The gauged Abelian symmetry does not represent the usual interaction between electric charge and the gauge boson. As a 
matter of fact, it is the additional global Abelian phase symmetry that happens to be related to the electric charge. 
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where A± are complex scalars, ip± are Dirac spinors and F± are auxiliary scalar fields that obey the following 
supersymmetric transformations: 

8A± = eV± (18) 
5ip± = sF± + ie^d^A± (19) 
5F± = iTfd^ . (20) 

Notice that, by introducing A± and ip± together in the same superfield, the scalars cannot be neutral 
under the global fermion number U(l) symmetry [341 136] mentioned in the previous Section. 

In (2+l)-D, N=1-SUSY accommodates neutral matter in real scalar superficlds. To describe charged 
fields, we have to combine two real supermultiplets. 

In the Wess-Zumino gauge, the gauge superconnection, T ai is written as 

T Q = i (j fl 8) a + 99X a and T a = -i (9^) a A^ + 96\ a , (21) 

where A^ is the gauge boson and A a is its partner, the gaugino (Majorana spinor). Defining the "field- 
strength" superfield W a as 

W a = -~D b D a T b , (22) 

with covariant derivatives given by 

D a = d a - i (>f6) a and D a = -d a -i (frf) a , (23) 

we obtain 

W a = \ a + TfthF^ - l -99^ ab (d^X b ) (24) 

and 

W a = A a - AKbF.u - l -ee (3 m a 6 ) , (25) 

where = \ [7^,7"] are the generators of the Lorentz group in (2+1) space-time dimensions. 

The gauge covariant derivatives that act on the matter fields with opposed L/(l)-charges, and 
are respectively given by 

V$± = (D a T iqT a ) $± and = (A, ± iqT a ) *L • ( 2 6) 

Using the definitions previously given for the superfields in Eqs. (fTT)) . (|2"TT) . (j2~4")l . (|2"5j) and the gauge 
covariant derivatives given in Eq. (|26l) . we can build the supersymmetric T3-QED action [45] : 

5 r3 _ Q ED = J d 3 xd 2 9 I —WW + (V$^)(V$+) + (V«l) (V*_) + 2m - | . (27) 

In terms of component fields, the action given by Eq. (|27l) may be written as 

SV 3 -QED = J cPx^iJ-fd M X-^F MV F' w + 

- A* + UA+ - A*_UA- + iij + ^d^ + + Up^d^- + F^F+ + F*F_+ 

- qAf, (i> + ^iP + - ^p_Y l ^P- + iA* + d^A+ - iA*_d^A_ - iA+d^Al + iA-d»A*_) + 



- iql A + tp + X - A^tp_X - A* + Xip+ + + q A^M 1 \A\A+ + A*_A_ j + 

-m(i} + tp + -i>_ip- + A* + F+ -A*_F_ + A + F$_ -A_F* ] 1 . (28) 



As one can easily read from Eq. (|28[) . the two-flavoured fermionic sector presents diagonal minimal 
coupling with the gauge boson, in correspondence with the first two terms displayed in Eq. (|15[) . Nevertheless, 
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the Yukawa-like interaction sector of Eq. (|15l) does not show up in the above T3-supersymmetric Lagrangian, 
even if we take the auxiliary fields F + and F_ on-shell. As a matter of fact, a $ 4 -sector has to be supplemented 
in order to provide the theory with the proper vertices, and such an extension must adopt the symmetry 
content of Lagrangian (f28f as a guideline for its construction. Thus, we re-state the invariance of the 
supersymmetric T3-Lagrangian with respect to an Abelian gauge transformation associated to the potential 
and to the field-strength F^ u , and we render evident a symmetry with respect to the following discrete 
parity variation: 

$± -> -$ T , (29) 

as space-time transforms according to x 1 , x 2 ) — > (x°, —x 1 ,! 2 ). 

The superfield variation expressed in Eq. (|29|) has its component-field counterpart spanned as 



V4 = WiM (30) 



A' ± — A T ; 
F± = Ft, 

where the 4-component spinor representation of parity reads 

*' - 7 « 2 )*- < 31 » 

with r Par defined in the context of a 4-spinor Pisarski representation, related to the Jackiw-Pi representation 
as drawn in the Appendix. 

To complete the variation rule for the supersymmetric action, we shall comment on the fact that the 
Grassmannian (2+l)D-intcgration measure has an odd character with respect to parity, 

d 2 e = dSde ^> d 2 e' = - d 2 e. (32) 

This stems from imposing homogeneity on the effect of parity action onto superfields $4-. The scalar 
superfields $± accomodate a contracted two-component fermionic term in their expansions, namely, 6ip±. 
As ip± transforms according to Eq. (|3"0")l . 9 is required to change tcH j°"f 2 9, up to a phase, as the term (9ip±)' 
must still be a scalar. So, regardless the phase, d 2 9 is necessarily odd under parity. As a consequence, the 
superpotential to be added to the theory has to be odd as well (note that the bare mass term proposed in 
Eq. (|27|) already shares such a property). 

4 The <I> 4 -sector 

The most general local U(l) and parity-invariant supersymmetric action with a $ 4 -sector reads: 

5<e,4 = fjd 3 xd 2 9 ($ 1 f $+) 2 - ($ f _$-)" 

+h J d 3 xd 2 9 ($ 1 f $_$ + $ + + $+$ t _* 1 f $ 1 f ) - (<^$+$_$_ + $_$^$t$t) (33) 
= Sf + Sh, 

where / and h are (real) coupling constants (to be associated later with Jackiw-Pi's g present in Eq. (fTS]) ). 
The quartic term with coupling constant h explicitly breaks the global U(l) fermion number symmetry. 
Nevertheless, we decide to keep this term in order to describe possible regimes where the chiral interaction 
dominates over the forces dictated by the global symmetry. The chiral symmetry has a dynamical character 



3 Such a variation also matches consistency as one considers both super-translations, — > + ie-y^O; 9 — > 8 + e, and parity, 
x' 1 — > —x 1 ; x°,x 2 left unchanged. Also, the one-flavor content of the Grassmannian sector of superspace, as one deals with an 
N=l theory, requires the variation of 9 to be associated to a 2x2 block of r Par . As a room for further extensions, we comment 
on the fact that the model here proposed is obtained after a truncation of an N=2,D=(2+1) descent of a N=l,D=(2+2) 
supersymmetric theory I45| . 
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in that it dictates a gauge interaction; on the other hand, the global fermion number appears to have only a 
kinematic character. The latter classifies the states and field configurations (e.g., vortices in Refs. [3"4l 155] ) 
without, however, introducing gauge-type interactions. Nothing prevents us from setting h — whenever we 
wish to recover results for which fermion number conservation is mandatory. 
The component-wise $ 4 -action expression follows: 



S f = d 3 xC f , 



with 



(F + A* + + F;A + ) \A+ 1 2 - (F_A*_ + F*_A_) \A_ | 2 + 2 (\A+\ 2 ^ + ~ |A_ | 2 ^_) + 



+- (Aii> + r+ + Afr + ^+ - A 2 _^_n 

Also, Sh = J d 3 xCh hosts the Lagrangian density 
-U/h = 



A*J^J-) 



(34) 



+ \A+\* [ A_F+ + A*_F; + -(A + F- + A* + F1) + ^+i>- + 1>-1>% ) + 



- \A- 1 2 + A* + F1 + - (A_F+ + A*_ F£) + + ^+4>- 

+ \ {Al4> + 4>- + A* + 2 ^_^ + ~ A 2 Jf_i> + - A* 2 ^-) + 
+ (A-A+ + A*_A* + ) (V>+V+ - + 
+ 1 [A* + A- - + A*_A+ - Wi'-)] + 

+i (A\A-F* + + A* 2 A*_F+ - A 2 _A+F*_ - A* 2 A* + F-) , 



(35) 



where the two-component charge- conjugated spinors are defined by ^£>S_ = ia x tjj^_ and = ia x ij)*_. 

As we aim to end up with an action on-shell for the auxiliary sector, the field equations for F± are 
obtained: 

SS 



SF+ 

SS 
6F1 

SS 
~SF* 

SS 

~sK 





— mA* + 


-fA\\A+\ 2 


- h 


\A+ 


2 A_- 


r \ (A* 2 A*_ 


-LA-| 2 ^-) 


= 


F_ 


+ 


+ fA_\A_\ 2 


+ h 


\A- 


2 A\- 


f \ (A 2 -A + 


-\A + \ 2 A%) 


= 


F+ 


— mA-\- 


-fA+\A+\ 2 


- h 


\A+ 


2 A*_- 


r \ {A\A- 


-\A-\ 2 A*_) 


= 


F* 


+ mA! 


+ fA*_\A-\ 2 


+ h 


\A- 


2 A + - 


f - (A* 2 A* + 


-\A + \ 2 A + ) 


= 



(36) 



The physical scalar sector of the complete action (that results from summing up Eq. (|27l) and Eq. 
reads as follows: 



V = F^F + + F1F_ = 

= m 2 (|A + | 2 + |A_| 2 ) + 2m/ (\A + f + |A_| 4 ) + mh (|A+| 2 + |A_| 2 ) (A+A_ + A* + A*_ 

+ ^ (|A+| 4 + |A_| 4 ) (A+A- + A* + A*_) - fh\A + \ 2 \A_\ 2 (A + A_ + A* + A*_) + 



f 2 (\A, 



| 6 + |A_| 6 



A A 



\A- 



' 2 )( 



4 I +11 



A, 



\A- 



h 2 . , 



AiAl + A\ 2 A* 2 



\A. 



+ 



(37) 
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where a minus sign has been omitted to allow for the direct identification of the potential V. 

At this point, we would like to reassess the procedure based on parity-symmetry we have followed to 
propose the <£> 4 -sector. As a matter of fact, the mass parameters present in extended models that follow 
original Jackiw-Pi's chiral gauge theory [35J [36] stand for v.e.v.'s of scalar fields, and a particular field is 
taken to be odd with respect to parity. Coherently, the corresponding mass parameter flips its sign as parity 
acts upon fields. If we had assumed the parameter m displayed in Eq. (f27|) to be odd as well, the parity 
transformation here proposed would not leave the action invariant, as it would also change by a sign. Instead, 
we take the viewpoint of considering the mass parameter as a fixed "bare" mass, here included for the sake 
of generality, left unchanged by parity as the physical mass would get contributions from non-trivial minima 
of the power-counting renormalisable ^-potential (mind Eq. (13710 . The corresponding critical values of 
proper combinations of A-fields should inherit the odd-parity property proposed by Jackiw et al. Also, as 
far as the parameters / and h are concerned, those are taken invariant under parity as well, for the premise 
of their identification with v.e.v's of scalar fields would lead to a non-renormalisable ^-potential, which we 
better avoid. 

For the simple case h = 0, the associated potential bears four minima: 



where the entries refer to values of \A+\ and \A_\, respectively. It is important to mention that, in order 
that spontaneous symmetry breaking takes place, it is necessary that / > and m < 0. Supersymmetry is 
preserved, since the potential in Eq. Q37p is zero for anyone of the minima above. 

If we consider the mass matrix for the whole set of original fermionic fields (A, tp + , ip'L, ipZ.) and 
evaluate its eigenvalues on, say, the asymmetric ()40[) minimum, we are left with the following results: two 

degenerate null outcomes; one — m eingenvalue; and two masses given by — y(l ± J~2 — ^j) [notice that, 

since q is real, m < and / > 0, these masses are real]. So, massive Dirac fermions appear in graphene (as 
in Ref. [33]) by means of the simplest self- interaction scenario. In addition, the corresponding eigenvectors 
are easily obtained and, by virtue of the breaking of gauge symmetry by the v.e.v. of the scalars, they 
obviously do not possess anymore a definite chiral charge Q Qualitative similar results are obtained when 
the two other non-trivial minima are considered. Massive charged fermions in graphene, after supersymmetry 
breaking takes place through some mechanism (we do not address to this issue here), might combine into 
scalar bound states and their attractive interaction may reveal non-trivial and interesting effects of the 
doubly-charged scalars. 



(0, 0), trivial case 



(38) 




(39) 



(40) 



(41) 



4 Of course, although not explicitly mentioned therein, the same occurs with the Fermi fields *S/± in Ref. 1341 whenever the 
scalar tp acquires a non-trivial v.e.v. 
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5 Conclusion and Perspectives 



Back to the component-wise expression for the spinor-scalar 3- and 4-vertex sectors of the complete super- 
symmetric action, we write 



d 3 x < -iql A+tp.X - A-ip_X - A* + Xip+ + A*_Xi/j- 



La% 2 + ±a%a 



2f\A+\ 2 + h(A + A_ + A* + A 
- 



2/|j4_ | + h (A + A- + A* + A 

+ h \a_\ 2 - - (a\ - a*J 2 ) - | (a; 2 - a 2 _) ^_v>+ 



La 2 4-4,4* 

2 + + 2 + - 



—A* 2 + - A* A< 

2 - 2 " H 



-A 2 + - A- A* 

2 2 4 



^1 - h\A + \ 2 (i)^ + i)_^l) 
h 



V+Vh- 

(42) 



The content of the above Lagrangian density together with the fermionic minimal couplings with the 
gauge boson displayed in Eq. (|28|l demonstrate that the supersymmetric T3-QED theory given in Eq. (|2T|) 
supplemented by a 4> 4 -term given in Eq. (|33|) provides a theoretical framework that extends Jackiw-Pi's 
original chiral gauge theory [34], represented here by Eq. (fT5|). Furthermore, room is also available to allow 
for accommodating extended models |35[ 36{ in a supersymmetric scenario, with the expectation that it 
may describe some of the physical features of graphene. Of course, some identifications are necessary: for 
example, in a crude comparison, one would identify Jackiw-Pi's scalar field ip with squares of the type A 2 ^, 
A 2 _ (and complex conjugates) which arise from our quartic superheld action, Eq. (|33[) . To exactly identify 
terms and degrees of freedom one should establish the proper combinations of fermions that diagonalize 
the mass matrix (eigenvectors) upon a particular choice of scalar fields configuration that minimizes the 
potential. This has been done at the end of the previous Section, for the case h = 0, which is the one to be 
considered for a comparison with Ref. [31], since, as already mentioned after Eq. (|33j) . the h-texms break the 
global ?7(1) symmetry, related to electric charge. All the results obtained are compatible with Ref. [54"] . 

We would like to stress — and this is a non-trivial remark — that the supersymmetric T3-QED related to 
the chiral gauge theory is nothing but a planar descent of an (N=l)-supersymmetric Abelian gauge theory in 
(2+2) space-time dimensions [15] . It remains to be understood — and, hopefully, we shall also be clarifying 
this issue soon — how to relate graphene physics to the dimensional reduction of a gauge theory set up in 
Atiyah-Ward space. This point sets out as a rich problem to be understood. 

A comment on the gaugino is in order. It is an additional field with respect to the work of Ref. [34] . 
To properly interpret the physical role of this (neutral) fermionic field, one has to consider that it mixes 
with the charged fermions whenever the complex scalars acquire non-trivial vacuum expectation values, as 
shown at the end of the previous Section. In this case, the A-field itself does not describe a mass eigenstate 
and, therefore, it cannot be interpreted as a physical excitation. So, after choosing a particular vacuum, 
it is possible to diagonalize the mass matrix in order to identify the correct combinations of the fermionic 
fields that correspond to the mass eigenvalues. This procedure is outlined in the previous Section, for the 
case h = 0. We could also check that SUSY is not broken (for h = 0). To the best of our knowledge, these 
properties cannot yet be accommodated into the experimental results of graphene. 

The physics of graphene has a very rich phenomenology. Supersymmetry may indeed open up a number 
of interesting issues in connection with the couplings between the fermionic and bosonic excitations that 
may arise in the description of this monolayer system. The spontaneous breaking of N=l-supersymmetry in 
(2+l)-D is less trivial than the same question in (3+l)-D, since N=l does not exhibit a complex structure 
underneath, as planar N=2-supersymmetry does. So, we conclude this work with the perspective to reassess 
the problem of spontaneous supersymmetry breaking in (2+l)-D by adopting graphene as a physical system 
to apply our results. In particular, supersymmetry breaking should help us to understand and classify the 
rich structure of couplings in graphene. 
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A Appendix 

A.l General Considerations 

All the gamma matrices respect the anticommutation and Hcrmitization relations respectively given by 

{7^ , l v } = 2rT, 7 Mt = 7V7 , (A.l) 
where if* = diag(l, — 1, —1). The other gamma matrices will always have the following properties: 

75 2 = 1, 75 f = 75 ; 7 3 ' = -1, 7 3t = -7 3 (A.2) 
and the additional anticommutation relations given by 

, 75} = , {7" , 7 3 } = , { 7 3 , 75} - . (A.3) 

A.2 Jackiw-Pi representation of gamma matrices 

The 4x4 gamma matrices in Jackiw-Pi rep. |34) . which will be called y\, can be expressed in terms of 
gamma matrices in minimal 2x2 representation, which we will represent by 7^, as 

u _ / 7^7° 



^ = { W ) (A - 4) 

with ji spanning only the values (0,1,2) and 

7 M = (a z ,ia y , -icr x ) . (A. 5) 

The additional matrices 7/ and 7j, are given by 

,J-( 1 °\ „,3_,.,o„i^_ f0 -7° 



7 5 J =^ 75 = nhWjrl - [ ; J. (a.6) 

Since 7j is not associated to the space-time, it will behave as a second 7 5 -type matrix, that is, it anticommutes 
with all the gamma matrices associated to the space-time. 

A.3 Pisarski representation of gamma matrices 

The 4x4 gamma matrices in Pisarski rep. |47j . which will be called 7p, can be expressed in terms of gamma 
matrices in minimal 2x2 representation, which we will represent by 7^, as 



-7^ 

with /J spanning only the values (0,1,2) and 

7 M = (a z ,ia y , -icr x ) . (A 
The additional matrices 7^ and jp, are given by 



Tf^lj ), 1% = i'fplWpli = { j oj- (A.9) 

Since jp is not associated to the space-time, it will behave as a second 7 5 -type matrix, that is, it anticommutes 
with all the gamma matrices associated to the space-time. 
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A. 4 The transformations that relate Jackiw-Pi. rep. to Pisarski rep. 



With the aid of the matrix S = I _ 0)1 we can write 



7£ = StftfS 
7 5 P = S^jS 
1p = Sj 3 jS . 



(A.10) 
(A.ll) 
(A.12) 



Eq. (|A.12|) follows from 7p and 7j definitions given, respectively, in Eqs. (|A.9[) and (IA.6I) and, also, 
from the transformations given in Eqs. (jA.101) and (| A. 1 1|) . In the sequel we will write down the inverse 
transformations, 
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